INTERSTELLAR HELIUM IN INTERPLANETARY SPACE 

William C. Feldman , /. /. Lange , and F. Scherb 

ABSTRACT The velocity distribution function of He* in the solar wind at 1 AU is calculated with the 
assumption that the source is photoionization of a cold (T = 100° K), neutral interstellar 
wind. If the spiral magnetic field is noise free, the velocity distribution is diffuse and 
would not produce a peak at 4(E/Q)jj in an E/Q particle spectrum. If the velocity of the 
interstellar wind with respect to the sun lies in the ecliptic, a large variation of the He* 
number density with respect to ecliptic longitude is expected. 


INTRODUCTION 

A peak appearing at four times the energy/ charge (E/Q) 
location of the hydrogen peak in two Vela 3 electro- 
static analyzer spectra of solar wind ions was interpreted 
as being mostly due to He* [Bame etai, 1968]. The 
He* abundance relative to He** was measured to be 
3X1 0 -3 ±50 percent (Bame, p. 548). Two other spectra 
published by Bame etai [1968] yielded only an upper 
limit of 10' 3 for n^ Q +(r e )]n}{ e ++(r e ), t ^ e ratio of num- 
ber densities at r = r e = 1 AU. A more recent attempt 
[Bame etai , 1970] to determine the mass and charge 
composition in the solar wind attributed a peak located 
at four times the E/Q location of H* as probably due to 
Si 7 * instead of He*. 

The Vela 3 observation cannot be explained simply 
since n^ e y(r e )ln^ + y(r e ) is expected to be of the order 
of 3X1 0^ 6 for a 10 6 °K corona [Tucker and Gould , 
1966]. Hundhausen etai [1968] attempted an expla- 
nation in terms of a model proposed earlier by Patterson 
etai [1963] . This model considers the solar wind pro- 
tons converted by charge exchange with galactic hydro- 
gen into hot, isotropic neutral hydrogen atoms at an 
assumed spherical shock boundary surrounding the 
heliosphere. These fast neutrals flow in toward the sun 
where a small fraction charge exchange with the solar 
wind He** to produce the He*. They find that if 
1 percent of the solar wind protons find their way back 
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and interact within 1 AU, enough He* could be 
produced to explain the Vela 3 data. Holzer and Ax ford 
[1970] tried to explain the ratio of 3X 10 -3 by assuming 
that the solar system is moving through an interstellar 
gas cloud. Neutral He would be gravitationally focused 
by the sun, photoionized, and then incorporated into 
the solar wind as He*. 

Rocket observations of diffuse extreme ultraviolet 
radiation (EUV) [Young etai, 1968, Johnson etai , 
197 la] from the night sky yielded an upper limit of 2 R 
of 584A radiation [Johnson etai , 1 97 lb] .If all diffuse 
584A radiation were to be interpreted as due to the solar 
584A line backscattered by neutral He in interplanetary 
space, then 2 R corresponds to ^pj e (°°) — 3X 10 -2 cm -3 . 
This is not a restriction to the Holzer and Axford 
[1970] interpretation of the Vela 3 measurements. 

In this paper we consider the orbits of interstellar 
He* ions injected into an ideal (noise-free) Archimedean 
spiral magnetic field. The He* ions originate from 
infalling interstellar wind neutral He atoms that are 
ionized in the interplanetary medium. We thus deter- 
mine the velocity distribution function of He* expected 
at 1 AU as a function of azimuthal position 0 in the 
ecliptic for different interstellar wind parameters. This 
study also gives orbits of atoms or molecules that escape 
from solar system bodies and are ionized in inter- 
planetary space. An example of such a source could be 
neutral hydrogen ionized beyond the disturbed interface 
between the Venusian ionosphere and the solar wind. 
Detection and identification of such particles could yield 
information concerning atmospheric abundances and 
thermalization processes in the solar wind. Finally, we 
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compare the various interpretations of infailing neutral 
He within the interplanetary medium. 

PITCH ANGLE SCATTERING BY MAGNETIC 
IRREGULARITIES IN THE SOLAR WIND 

Before proceeding to a discussion of the orbits of very 
low-energy charged particles injected into the solar wind, 
we consider the magnitude of the perturbation caused 
by cyclotron-resonant pitch-angle scattering from mag- 
netic irregularities convected by the solar wind. The 
problem with relevance to the low-energy particles of 
interest here has been treated theoretically by Sturrock 
[1966] and Hall and Sturrock [1967] , and applied by 
Barnes [1970] to fast particles accelerated at the earth’s 
bow shock and moving upstream in the solar wind. 
Assuming the dominant scattering process is pure mag- 
netic cyclotron-resonant pitch-angle scattering, the par- 
ticle velocity pitch-angle distribution function F is well 
approximated by a diffusion equation of the form 

f 1 ■ a * v u iz) f a s? & h eG f ] 

where Fy is the particle’s velocity component parallel to 
the unperturbed magnetic field f? a , and 
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where v = eBjlnmc is the cyclotron frequency, and 
P(y) is the power density spectrum of the transverse 
magnetic field fluctuations evaluated at the frequency v\ 
which is Doppler shifted into the cyclotron frequency in 
the rest frame of the particle. To obtain an order of 
magnitude of the time scale r for the disruption of a 
beam of He + particles moving relative to the solar wind, 
we approximate G as a constant and set* = cos 6 so that 



In the frame of reference moving with the particle 
beam 

CO 

F(x,t) exp [-1(1 + l)GtPi(x)} 

l=o 

where the Pj(x) are /th-order Legendre polynomials. 
Thus the decay time r/ is given by 

_ 1 
T/ /(/ + l)G 


Power spectra of the interplanetary magnetic field 
have been reported by Coleman [1966] for the 
Mariner 2 magnetometer, by SiscoeetaL [1968] for the 
Mariner 4 magnetometer, and by Sari and Ness [1969] 
for the Pioneer 6 magnetometer. Comparison of these 
results indicates that the power level at a given fre- 
quency can vary significantly depending on solar 
activity. Since the He + particles to be discussed later are 
for the most part moving away from the sun, a conser- 
vative choice of frequency for evaluating the power 
spectrum is the cyclotron frequency which for He + in a 
6y magnetic field (at 1 AU) is 2.3X1 0~ 2 Hz. Choosing 
P(2X10“ 2 Hz) = SXKTV/Hz from Mariner 4 data for 
quiet periods yields 

2G = 2XKF 6 see 1 

and consequently Ti= 5X10 5 sec. But 5X 10 s sec is 
about equal to the expected ion transit times from 
r = 0.1 AU to 1 AU if the He radial velocity at the point 
of ionization is away from the sun. 

This estimate for r is probably uncertain by at least 
two orders of magnitude due primarily to the approxi- 
mate nature of the diffusion equation, the evaluation of 
P at 1 AU and the choice of / = 1. Of course, more 
peaked pitch-angle distribution functions are character- 
ized by higher l values and so are damped more quickly. 
Most important, because it is unknown experimentally, 
is the radial dependence of the power spectrum evalu- 
ated at the Doppler shifted gyrofrequency. Evidence 
from Mariner 4 magnetic field data [Coleman etal , 
1969] indicates that the power density for a given ‘fre- 
quency generally decreases roughly as r~ 2 for 
1 AU <r< 1.43 AU where r is the heliocentric distance, 
but as r increases, B decreases and so does the gyro- 
frequency. The radial dependence relevant to our 
problem is that of the power density spectrum evaluated 
at the Doppler shifted gyrofrequency. From figure 14 of 
Coleman etal [1969] it is seen that this is almost 
independent of r for 1 AU O < 1 .43 AU. In any case, it 
is expected to be less dependent upon r than the power 
density spectrum evaluated at a fixed frequency. 

Considering these uncertainties coupled with the 
measured variability of the solar wind, one can 
conclude neither that significant pitch-angle scattering 
occurs in transit to 1 AU nor that the He + orbits remain 
unperturbed by magnetic field fluctuations. Further- 
more, as pointed out by Barnes [1970] , refers only to 
the disruption of the beam by hydromagnetic waves but 
not to thermalization. The time scale for thermalization 
is expected to be longer. 

It is thus useful to determine the unperturbed orbits 
in an ideal spiral magnetic field as a function of injection 
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position and velocity. It is expected that particles 
injected at small r will undergo some degree of 
pitch-angle scattering, whereas those injected around 
1 AU will follow unperturbed orbits. One can then com- 
pare future measurements of the He* velocity distribu- 
tion function with both the thermalized and unper- 
turbed limits to obtain information concerning 
interstellar wind parameters as well as wave-particle 
interactions in the solar wind. 

UNPERTURBED ORBITS OF A TEST PARTICLE 
INJECTED INTO THE SOLAR WIND 

The coordinate system relevant to this discussion is 
shown in figure 1 . The sun is assumed to move with 
velocity (Voo) relative to a cold ( T= 100° K) interstellar 
cloud such that in the frame of reference with the sun at 
rest at the origin, neutral He atoms appear to be moving 
parallel to the YZ plane at an angle 0 with respect to the 
Z axis. The positive Z axis is oriented northward perpen- 
dicular to the ecliptic. The XY plane is coincident with 
the ecliptic. Infalling He atoms follow hyperbolic trajec- 
tories before being ionized by solar EUV. After ioniza- 
tion they are accelerated by a radially convected spiral 
magnetic field. Although the spiral magnetic field 
pattern is cylindrically symmetric, the subsequent 
charged particle trajectories are not, due to the 
asymmetric initial injection velocities. 



Y ECLIPTIC PLANE 


Figure 1 Interstellar infall coordinate system with 
(V^ in the ecliptic. 


Each He atom hyperbolic orbit is planar. Every point 
in the orbit plane can be specified equivalently by r . the 
radial distance from the sun, and 5 , the angle between r 
and <Voo>, or by two impact parameters p t and p 2 identi- 
fying both a direct and indirect orbit, respectively, 
which pass through that point. Fahr [1968] obtained an 
expression for the neutral number density n^(t) at any 
point for T- 0° K using the continuity equation. 
Rewriting his results we obtain 


«£j e ( r ) = «i exp 




«2 " «i - « He (°°) 

2 1 1/2 

rsinfi /rsin§\ 

Pi = —^~ + — ) + rb(l + cosS) 

p 2 = p y - r sin 5 
GM 

b= '2 

(vj 2 

and w j| e (°°) is the number density of neutral He at r = *>, 
P(r e ) is the probability per sec for ionization at 
r e ~ 1 AU, G is the gravitational constant, and M is the 
solar mass;P(r e ) is variable, but is chosen here as 

P(r e ) = 6.8X1 0“ 8 sec^ 

In figures 2 and 3, we plot two cuts through the nor- 
malized neutral He number density distribution 
The radial dependence of 
(r)/«„ C 010 ) for (V^ = 20 km/ sec oriented in the 
ecliptic (0 = 90°) is plotted in figure 2 for ecliptic posi- 
tions oriented at right angles to <Voo>- One sees that, 
beyond 1 AU, ^jjgW/^jjgC 00 ) * s essentially constant. On 

the other hand, the azimuthal dependence is dominated 
by a dense tail in the antiapex direction. This feature is 
characteristic of the gravitational focusing problem; it is 
real and limited only by the temperature of the infalling 
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Figure 2 Normalized neutral He number density as a 
function of heliocentric distance: = 20 km /sec, 

®=90°, 0- 0° or 180°. 



<3> (DEGREES) 



Figure 4 Velocity space coordinate system relevant to 
calculation of the neutral He number density for 8 = 0. 



Then, 


Figure 3 Normalized neutral He number density as a n ne <r,o)‘ 
function of azimuth: (V^ = 20 km/sec , T= 100 K, 
r = 0.5 A U. The He density plotted in the neighborhood 
of 0- 270° is the sum of contributions from both direct 
and indirect orbits . 


n°^(r,o)F(y o ^V < J dV^ sin 6 dS dr] 


gas. For 8 = 0 and 7= 0 K,p« ~p 7 ; hence, both a j and ^ A 0 / <>\ * + 

^ , To find n" (r, 5), note that 

n 2 .~ °°. For T different from zero but such that ™ 

m<Koo> 2 l2kT » 1 , we can obtain an analytic expression 

for n , and n 2 at 5 = 0 by integrating n°_ (r) weighted by lim p 1 = ^ 2br) lim (pi ~ = 

an isotropic maxwellian velocity distribution function 

over V^. The velocity-space coordinate system for this hence 

calculation is outlined in figure 4. 

For 5 « 1, we can write a maxwellian distribution ( r g\_ 2. n / I - 7 ir e 2 P(r e ) j 

in the form He ’ S He \ry I (2br) in <f^o> 
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and 

"He^- 0 ) 



Here we have extended the limits of integration for 5 to 
+°° since mV 0 o(B 00 )l2kT» 1 and the exponential will 
cut off the integral. The 5 integral can then be evaluated 
exactly. Integration over is straightforward. The 
result is 


Ds MVJ 2 

21kT 

The neutral He-normalized number density distribution 
is plotted in figure 3 as a function of ecliptic azimuth 0 
for (Voo) = 20 km/sec, 100° K, and r = 0.5 AU. The 
upper curve refers to direct orbits (characterized by 
impact parameter p t ), and the lower curve refers to 
indirect orbits (characterized by impact parameter p 2 ). 
The He density plotted in the neighborhood of 0- 270° 
is the sum of contributions from both direct and indirect 
orbits. Two important characteristics of the curves 
should be noted. First, in the forward hemisphere direct 
orbit densities dominate indirect orbit densities due to 
the much longer angular distance traveled by the latter. 
Second, there is a marked anisotropy (here 67:1) of 
number density in backward and forward directions 
resulting from an enhanced tail with FWHM = 3.1° for 

r= ioo° k. 

We now specify the ion position in the usual ecliptic 
polar coordinates r, 0, 0 where r is the radial distance 
from the sun, 6 is the polar angle with respect to the 
north ecliptic, and 0 is the azimuthal angle with respect 
to the x axis. We further define X as the angle between r 
and the magnetic field B. Then the Archimedean spiral 
magnetic field is given by 


, , , JnGMm\ in [-tt r e 2 P{r e )~\ 

"H *"H«M Uf) 




Off-axis values of «jj e (r,8) were constructed using the 
parabolic-interpolation formula 

where D is determined by the imposed constraint 

B Hefc°Xl - D8 C 2 )- j- K H e(°°) ( 7 ) X 

r-nr e 2 P(r e y] 

at some angle 8 C . This yields 



and 


b g = o 

2 } (2) 
B <P=- B a(j) tan X 

2 

B = B a( j)/ cosX 

where cos X = [1 + (rco sin OjV sw f ]~ 1/2 , the angular 
speed of the sun go = 2.9X10" 6 rad/ sec, B a is the radial 
component of the magnetic field at some reference 
radius#, and ^sw is the radial solar wind velocity. 

The equations governing the motion of a charged test 
particle in a slowly varying electric and magnetic field in 
the guiding center approximation were developed by 
Northrop [1963] . The motion can be separated into 
three components: (l)’the superposition Vp of various 
drift velocities perpendicular to the magnetic field, (2) a 
circulating velocity Vi about the magnetic field, (3) a 
velocity Fjj parallel to the magnetic field. We now treat 
each of these components in turn specializing to the 
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above radially convected spiral magnetic field pattern. 

The only perpendicular drift velocity of consequence 
is that due to the appearance of an electric field in the 
frame of reference at rest with respect to the sun. Here 
Vf) = V sw sin X with direction perpendicular to B but 
confined to the cone specified by 6 - constant. 

The component of motion circulating about B is 
determined in the adiabatic limit (which fits our 
problem) by the conservation of the magnetic moment 
/ u, given by 



(3) 


Combining equations (5) and (6) to eliminate the time, 
we obtain (after some algebra) an equation for Fy in 
terms of r alone: 



= dr [£ T^ (1 +COs2 X) - 7? + (F sw 2 ) 

Equation (5) can be integrated to yield 


sin 2 Xcos 2 X 


(7) 


The equation governing the motion parallel to the 
field direction given by Northmp [1963] is 


V r = (V || cos X + V sw sin 2 X)=± 


Vswffi) 

(i +n i/2 


mdV || 9 B 

— — - mg\\ + qE\\~n — +mV D = ( . F|| sin x + Fjw sin x C os X) 


9 e U 

9e u 


r .... 

11 

jr + 

V\\—+(V D -V>\\ 

(4) 

=(i +f 2 ) i / 2 L (i+f } T/(r) J 


where m is the mass of the particle, the magnetic field 
direction is specified by the unit vector ej|, gy is the 
component of the sun’s gravitational acceleration along 
ey, q is the particle charge, E y is the component of elec- 
tric field parallel to ey, and s is arc length along ey. In 
our problem i?y = 0 and 3ey/3£ = 0. Setting 3/3 s = ey • V 
in equation (4) and carrying through the algebra, we 
obtain 

dV \\ n B cos X - „ GM cos X 

~=- (1 + cos 2 X )-— -3 


F, w sin 2 X 

+ — [V sw cos 3 X- F||(l + cos 2 X)] (5) 


The total noncirculating velocity is then given by 


where 



Here f -cor sin d/V sw is a dimensionless radial distance, 
and all quantities appearing with a subscript zero refer to 
the initial conditions. The equation for the trajectory 
can be written in dimensionless form as 


yielding 


V^e+Vtf 


V\\ cos X + V sw sin 2 X 


df _fsin0F r sin 6 /(f) 

di> = ~$ =± [o+r 2 ) 1 / 2 */(f)] (11) 

where the minus and plus signs refer to motion in, 
toward or away from the sun, respectively. 

To complete the problem the initial velocity condi- 
tions at the point of injection need be specified. We 
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consider first the trajectories resulting from a zero initial 
velocity and then proceed to the interstellar wind prob- 
lem. 

Zero-Velocity Initial Conditions 

Here Fj| Q=0 and Fj_q = Vq = V sw sin X 0 since the 
orbits are initially cycloidal. Equation (1 1) was solved 
on the University of Wisconsin’s Univac 1 108 computer 
with a “Math-Pack” subroutine utilizing Hamming’s 
method of solving ordinary first-order differential equa- 
tions. The initial noncirculating velocity is perpendicular 
to B but away from the sun. Since B decreases as r 
increases, the conservation of the magnetic moment con- 
verts part of the circulating velocity into the parallel 
component. Specializing to the problem of injection into 
the unperturbed solar wind at Venus, the azimuthal 
angle as a function of radius characterizing the trajectory 
is plotted in figure 5. Components of the velocity as a 
function of r are plotted in figure 6. Again, 
V sw = 400 km/sec. Since both the spiral field pattern 
and the initial injection conditions are azimuthally 
symmetric about the sun, there is no focusing of the 
particle beam in the ecliptic. Thus it’s observed angular 
diameter should remain constant. Any deviation would 
then be due to interaction with magnetic irregularities 
convected by the wind. 



Figure 5 Azimuthal position of a charged particle 
relative to the azimuthal position of injection into the 
solar wind at Venus. 


Interstellar Infall Initial Conditions 

The neutral particle velocity at the point of ionization is 
determined by conservation of energy and angular 
momentum. For the indirect orbits 
v r =M1 + 2 b\r Q - (p 2 lr Q f ] 1/2 and V y = -VooP 2 lr 0 , 
where b-GM/V 0 2 . For the direct orbits 



Figure 6 Components of the velocity of a charged 
particle injected into the solar wind at Venus as a 
function of heliocentric distance. 


V r = ±Ko[l + 2 b/r - (pilr) 2 ] 1 12 where the ± refers to 
7 ^ Tmin = tan_1 (rPilb ) and Vy“ V OQ p 1 Jr 0 . ] The initial 
neutral velocity is Vyy = V r e r + V^ey yielding 
(F||)oo = Vyy ^ e|| and (V])oo = V/y- V N • ej| — V D . The 
effect of the asymmetry of these initial conditions with 
respect to azimuthal injection angle 0 O on the subse- 
quent trajectories is greater as r 0 decreases. This asym- 
metry is illustrated in figures 7 and 8 for injection from 
direct and indirect neutral orbits, respectively, at 
r 0 = 0.15 AU withVoo and r 0 in the ecliptic. The effect is 
most noticeable for the direct orbits characterized by 
90° < 0 O < 270°. As 0 O increases from 90°, the mag- 
netic moment decreases due to a cancellation effect 
reducing the magnetic mirror repulsion of the diverging 
field lines. Consequently, there is a region of impact 
angles around 0 O - 180° for which ionized He atoms 
have access to the photosphere. This assumes a noise-free 
Archimedean spiral model extending all the way in to 
the photosphere which is not the case in reality. Pro- 
ceeding further around in azimuthal angle, V r decreases 
to zero as well, yielding trajectories that circulate about 
the sun through large angles before making their way out 
to large r. One such trajectory not plotted in figure 7 
circled about the sun three times before reaching large r. 
Of course, such an orbit is unphysical due to interaction 
with magnetic irregularities but it illustrates the azi- 
muthal focusing characteristics of the problem. The 
indirect orbits are illustrated in figure 8. These are much 
more regular but still demonstrate a degree of azimuthal 
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Figure 7 He + orbits in an ideal spiral magnetic field 
following injection from direct neutral helium orbits: 
(V^) = 20 km/sec, 0 -90°, V sw ~ 400 km [sec , and 
r o = 0.15AU ; 



Figure 8 He + orbits in an ideal spiral magnetic field 
following injection from indirect neutral helium orbits: 
(Fqq) " 20 km/sec, 0 = 90°, V $w - 400 km/sec, and 
r Q = 0.15 AU 

focusing. The effects of the asymmetric initial condi- 
tions on the trajectories become less pronounced as the 
injection radius increases since the infall velocity 
becomes small compared to the drift velocity 
Yd = V sw sin X. 


We are interested in calculating the He + velocity dis- 
tribution function, which would be measured by a space- 
craft observing at 1 AU in the ecliptic. Toward this end 
we need to know the locus of all source points for par- 
ticle trajectories intersecting a given angular sector at 
1 AU. This was accomplished numerically for 20° wide 
angular sectors centered on 0 = 0° , 90° , 180°, and 270°, 
for (Voo) = 10, 20, and 40 km/ sec, and for angles of 
incidence 0 with respect to the north ecliptic of 0°, 45°, 
and 90°. All orbits penetrating to r Q less than 0.05AU 
were excluded from consideration. Typical configura- 
tions are illustrated for direct and indirect orbits in 
figures 9 and 10 for <Kx>)= 20 km/sec and 0 = 45°, As 



Figure 9 Locus of source points of He + for detection 
at 1 AU following injection from direct neutral helium 
orbits: (V^ = 20 km/sec , 0 = 45°, V sw = 400 km/sec. 


can be seen, the focusing effect of the interplanetary 
fields becomes more prominent at small injection radii. 
For example, the locus of points intersecting the 0 = 0° 
sector from r Q = 0.15 AU splits into two angular sectors 
at 0.15 AU due to the partial cancellation of the mag- 
netic moment coupled with the changing sign (from 
minus to plus as <j> 0 increases past the angle of peri- 
helion) of the initial radial velocity . 

To construct the He + velocity distribution function 
we must sum separately the contributions enclosed 
within each sectored locus diagram for both direct and 
indirect orbits. Denoting direct neutral orbit initial con- 
ditions by J - 1 and indirect neutral orbits by J = 2, we 
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t<v„> 


Figure 10 Locus of source points of Lie* for detection 
at 1 AU following injection from indirect neutral helium 
orbits: {VJ} = 20 km/sec, S = 45 ° , V gw = 400 krrjjsec. 


obtain from flux conservation 


where 


" 2 - = / 


r e dN + 


dN + n He( r O’ < l > o) p ( r e) (d<$>\ 

dr ° vy\r e ,<t> e ;r 0 ,d> 0 ) WJ 


( 12 ) 


(13) 


and n'J^(r e 4> e ) is the He + number density at the point 

(r e , <t> e )\ n^(r 0 ,<P 0 ) is the neutral He number density at 

the point (r 0 ,<j> 0 ); vf/)(r e ,<p e ;r o ,0 o ) is the radial com- 
ponent of He + velocity at (r e ,<l> € ) for those particles 
originating at(r 0 ,$ 0 ); and d<j> 0 !d<p e denotes the ratio of 
the azimuthal angular sector d(p 0 at (^ O ,0 O ) to the given 
sector d(j> e at ( r e 4> e ) into which d4> 0 is mapped. The 
integral is evaluated along the locus sector. 

The integrand of equation (12) is plotted as a func- 
tion of r Q and <j> e in figures 1 1 through 14 for (VoJ = 10, 
20, and 40 km/sec and 0 = 45° and 0°. No plots for 
0 = 90° are given due to the extreme anisotropy 
resulting from the presence of the dense, well-defined 
neutral-particle tail mentioned earlier. The definition of 
this tail depends strongly on the neutral gas temperature. 

For T= 100° K its contribution to depends 

critically on © since the FWHM of the tail is about 3°, 



r. (AU) 

Figure 1 1 Differential He + number density at 1 AU as 
a function of injection radius: 10 km j sec, 

0 - 45°; <j> e denotes the azimuthal position of 
observation in the ecliptic . 



Figure 12 Differential He + number density at 1 AU as 
a function of injection radius: (V^ = 20 km/sec, 
0 = 45°; <p e denotes the azimuthal position of 
observation in the ecliptic. 


Furthermore, a crude (for 3° FWHM variation) 15° grid 
was used in our numerical analysis to calculate 

n^ + (r e ,0 e ) and thus is not suited for calculation in this 
case. 

For 0 = 90° and T~ 100° K, one might expect to 
(J) 

observe anisotropies in «„ e +(^,0 e ) of the order of 50 
(fig. 3). 
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r* (AU) 


Figure 13 Differential He + number density at 1 AU as 
a function of injection radius: (V^) ~ 40 km /sec, 
&~45° ; <p e denotes the azimuthal position of 
observation in the ecliptic. 


Table 1 Integrated normalized He + number density 
n H ^ (r e^e)l n n Q (^) at 1 AU 


© = 45° 

(KJ 


n ,(r p ,<t> p )ln (°°) 
He + ee He 


(km/sec) 

tp e = 90° 

<p e = 180° (j> e = 270° 

<p e = 360° 

10 

3.8XKT 2 

4.2X1 O' 2 7.9X1 CT 2 

7.5XKT 2 

20 

2.7X1 (T 2 

2.8X10’ 2 4 .4X1 O’ 2 

4.1X1CT 2 

40 

2.4X10" 2 

2.4X1CT 2 3.0X10T 2 

2.5X1CT 2 

© = 0° 


<Koo> 

(km/sec) 




10 

5.0X1 cr 2 



20 

3.2X1 cr 2 



40 

2.5X1 O' 2 



<V®> (km /sec) 



r 0 (AU) 


Figure 14 Differential He + number density at 1 AU as 
a function of injection radius for 0 = 0°. 

The He + number density at 1 AU integrated over 

velocity space and normalized to/ijj e (oo) is summarized 
in table 1 for 0 =45° and 0° and<K»> = 10, 20, and 

40km/sec. It can be seen that increases 

with decreasing < Voo ) and that the anisotropy with 
respect to azimuthal phase § e at 1 AU increases with 


increasing 0 as well as decreasing <K»). These effects 
primarily reflect the behavior of the neutral He number 
density in interplanetary space, but are also affected 
somewhat by the greater anisotropy of initial velocity 
conditions as 0 increases to 90° . 

From equations (3), (8), and (9), we know the veloc- 
ity of particles at 1 AU as a function of their injection 
position. Combining this information with dN + /dr Q 
from equation (13) we can construct the unperturbed 
velocity distribution function. Figure 15 shows a projec- 
tion into the ecliptic of the position in velocity space 
occupied by ions with injection radii ranging between 
r 0 = 0.15 AU and 1 .0 AU. The actual positions in veloc- 
ity space are circles about B with diameters indicated in 
the figure. A spherical velocity contour of 20 km/sec 
about the solar wind velocity V sw = 400 km/sec is 
sketched in for reference purposes. The minimum veloc- 
ity angle with respect to the radial direction is 21° for 
the conditions (Voo) = 20 km/sec, 0 = 45°, and 
0 e =18O°. This minimum velocity angle for different 
sun-cloud configurations is tabulated in table 2. It is seen 
that this angle varies only slightly for markedly different 
incident conditions. The cone of occupation can change 
shape by at most 1 5 percent to 20 percent for the dif- 
ferent incident conditions considered here. But since the 
distribution function is so diffuse to begin with, this is 
not a major effect. 
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Figure 15 Projection of the Re* velocity distribution 
function into the ecliptic. Each line drawn perpendicular 
to B represents the diameter of a circle about B corre- 
sponding to He* ions injected into the solar wind at 
r- r Q . The circle about a radial velocity of 400 km/sec 
represents the projection of a sphere with 20 km j sec 
radius. 

DISCUSSION 

We proceed now to the question of whether the model 
of cold interstellar He infalling into interplanetary space 
is consistent with the Vela 3 observations of He*. The 
two relevant Vela 3 observational results are: (1) Most of 
the time, the ratio n^+(r e )jn^+i{r e ) isle$s t han'lCT 3 ; 
and (2) on rare occasions, n \{Q J r( r e)l n l\Q++( r e) * s 
observed to be 3X10“ 3 ± 50 percent [Bame et al, 
1968]. 

Furthermore, we require that the model be consistent 
with recent observations of extraterrestrial diffuse 
Lyman a radiation observed with instruments on OGO 5 
[ Thomas and Krassa, 1971; Bertaux and Blamont, 
1971] . To interpret these observations, Thomas [1971] 
considered two models of interstellar neutral hydro- 
gen streaming relative to the sun. His first model 


Table 2 Minimum velocity angle (3 of the He* velocity 
distribution function with respect to the radial direction 


0 = 45° 


(km/sec) 

0 e = 9O° 

0 e = 180° 

<t> e = 270° 

0 e = 360 c 

10 

20.0 

21.0 

22.1 

19.9 

20 

19.8 

20.8 

22.1 

19.7 

40 

19.4 

20.7 

22.0 

19.3 


0 = 0 ° 

(VJ p 

(km/sec) (deg) 

10 21.2 

20 21.1 

40 20.9 


required a cold (T - 100° K) streaming flow with 
n H( 00 ) = ^-^ cm_3 9 10 km/sec oriented in the 

ecliptic, and an isotropic galactic Lyman a background 
of about 2007?. This was considered unlikely since cur- 
rent theories of the interstellar medium usually associate 
high temperatures ( T ^ 10 4 °K) with a hydrogen num- 
ber density as low as 0.06 cnf 3 . Furthermore, a calcula- 
tion by Adams [1971] of the galactic Lyman a back- 
ground radiation yielded an upper limit of 200/? with- 
out including the effects of Lyman a absorption due to 
interstellar dust grains. Thomas’ second model requires a 
hot (T^5X10 3 °K) flow with «H(^)^0.12cm~“ 3 , 
< Foo> — 6 km/sec oriented approximately in the ecliptic, 
and no galactic Lyman a background. (The upper limit 
of 2 R of 584A radiation set by Johnson et al [1971ft] 
does not contradict the presence of neutral He 
consistent with the above models.) 

The model of infalling interstellar gas considered in 
the present paper is valid for a finite gas temperature T 
at r=°° as long as m(Voo) 2 12 kT» 1. For the cold 
model, m<Foo> 2 /2kr= 243 » 1 , and our results are 
applicable. But, for the hot model (K 00 = 6 km/ sec, 
r= 5000° K\m(Voo) 2 / 2 kT~ 1.72 and our expressions 
for the neutral He number density in the antiapex region 
break down. In the following discussion we therefore 
concentrate on the cold interstellar gas model; the hot 
model will be discussed in a later paper. 

As mentioned earlier, placement of (VJ in the eclip- 
tic results in large variations in n^+(r e ,<p e ). Figure 9 
shows that enhanced fluxes are expected for 
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270° < <t> e < 360°. 

In the neighborhood of 0 e = 270° the major contribu- 
tion comes from neutral He atoms ionized in the dense 
tail close to 1 AU (see figs. 3 and 9). These ions occupy 
the region at the base of the triangular cross section of 
the velocity distribution as indicated in figure 15. Since 
these ions travel a relatively small distance in the inter- 
planetary field before reaching 1 AU, one does not 
expect much pitch angle scattering. Hence these ions 
would probably not contribute to a peak in an EjQ par- 
ticle spectrum at 4 (E/Q)jj. 

In the neighborhood of <j) e = 0°, the major contribu- 
tion comes from neutral He atoms ionized in the dense 
tail close to the sun. These ions normally occupy the 
region of velocity space at the apex of the triangular 
distribution which is closest to the solar wind velocity. 
In view of the evidence for enhanced noise in the solar 
wind near the sun, a high / value characterizing the initial 
pitch angle distribution of He + ions, and the long travel 
times to 1 AU, it is expected that a good fraction of 
these ions could be incorporated into the solar wind to 
contribute to a peak at EjQ = 4 (E/Q)fj. 

One might also expect a contribution to a peak at 
4 (E/Q)ff for azimuths in the range 100° <0^ 240° 
(for radial ionization distances 0 O 0 < 0.2 AU). Here 
the magnetic moment is small so that He + ions penetrate 
deeply into the corona (see fig. 7). Such ions would be 
expected to be thermalized and therefore incorporated 
into the solar wind. 

Last, a thermalized He + contribution is expected at 
<j) e = 270° due to neutral He in the tail ionized at 
^<0.1 AU. This He + number density is given by 


270 “)- J 


0.1 AU 


P(r e ) 

n He( r O’ 270°) — dr Q 

Vsw 


so that we obtain 


« He+ (^ 27 ° 0) 


He 




5X10' 


1-4 


,++ve 


on the basis of the solar EUV data of Watanabe and 
Hinteregger [1962]. A more recent measurement by 
Hall and Hinteregger [1970] revealed that the original 
measurements were systematically high. Use of the new 
data would decrease P(r e ) and hence increase the neutral 
He density close to the sun. This would tend to make 
the He + density at 1 AU larger than the values calculated 
above. On the other hand, we did not include collisional 
ionization in calculating R(r e ). This becomes important 
close to the sun which would tend to decrease 

n He+( r e> 2 7Q°y 

Hundhausen et al [1968] proposed, as an explana- 
tion of the Vela 3 measurements, that backstreaming 
hot neutral hydrogen atoms produced by charge 
exchange of solar wind protons with interstellar hydro- 
gen at the heliosphere boundary could charge exchange 
with solar wind He ++ , producing He + . This mechanism 
requires the flux of incoming hot (<F,o>~ 400 km/ sec) 
neutrals at 1 AU (which become ionized before leaving 
1 AU out the reverse side) to be 1 percent of the solar 
wind flux of H + . This fraction F(r e ) can be estimated if 
we assume that all solar wind protons are converted to 
an isotropic flux of hot neutrals at the heliosphere 
boundary. Then if /(H + ) is the flux density of solar wind 
protons at 1 AU, the total flux of fast H atoms back- 
streaming through any radius is 4 nr e 2 f. The fast H atoms 
move along nearly rectilinear paths. Specifying the 
impact parameter of a backstreaming H atom by p and 
defining 5 as previously, the probability a of an interac- 
tion within 1 AU is 



Then 


0 

where R is the radius of the heliosphere. Integrating 
numerically, we obtain 


for «He(°°) = 3X1 0“ 3 cm -3 , <F<x>> =10 km/ sec, T = 
100°K, n He ++( r e) = 0.25 cm" 3 , V sw = 400 km/ sec and 
P{re)~ 6.8X10~ 8 sec -1 . This ratio is not inconsistent 
with the Vela 3 result, but Vela 3 was at 0 e = 210° when 
the He + was observed. 

One further comment with respect to the above cal- 
culations is in order. P(r e ) = 6.8X 10~~ 8 sec -1 was chosen 


F(r e ) = 0.088 (j) 

But R is expected to be at least 10 AU so that we obtain 
F(r e )<^ 10~ 3 . This is an order of magnitude smaller than 
required by Hundhausen etal [1968] to account for 
the Vela 3 measurement. 
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CONCLUSIONS 

We conclude that although the possibility exists that 
«He+( r e^e)/ w He ++ ( r e) in the solar wind is substantially 
higher than predicted from ionization equilibrium in a 
10 6 °K expanding corona, it could be highly variable 
and depend strongly on 0 e , the azimuthal position of 
observation about the sun. We have shown that if the 
spiral magnetic field is noise free, the velocity distribu- 
tion of He + at 1 AU is diffuse. Such a spectrum would 
not produce a peak at 4(ElQ)fi in an E/Q spectrum. If 
(Vko) is in the ecliptic, large variations in with 

respect to <j> e are expected. Since solar wind magnetic 
irregularities are expected to be enhanced close to the 
sun, certain <j> e values will be more favorable than others 
for observing He + in a peak at 4(67<2)jj in an EIQ spec- 
trum. For example, azimuthal positions in the neighbor- 
hood of 100° <240° corresponding to injection 

radii 0 < r Q ^ 0.2 AU receive contributions due to a can- 
cellation of the initial magnetic moment, resulting in 
deeply penetrating He + orbits. This may result in effec- 
tive thermalization. Another likely range of locations for 
observing thermalized He + ions is 270° <j> e 360°. 

A value of n ^ e +( r e)/ w ji e ++(^) in the solar wind of 
^5X1 0“ 4 can be predicted by the presence of a dense 
neutral tail due to interstellar infall. More He + is avail- 
able from the interstellar infall, but a calculation of its 
relative magnitude appearing in a spectral peak requires a 
knowledge of the radial dependence of the power 
density spectrum of magnetic irregularities in the solar 
wind at times just preceding the observation. 

It is interesting to note that the Vela 3 spectra with 
the anomalous peak attributed to He + were obtained on 
October 8. October 8 corresponds to 0 e = 21O° for the 
earth. This is in the neighborhood where magnetic 
moment cancellation is expected to produce deep radial 
penetration and hence thermalization. Of course, deep 
penetration could also result in a second ionization 
yielding He ++ . We feel that one should not place too 
much emphasis on a comparison of a single measurement 
with the ideas presented above. 
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